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A NEW METHOD IN BOUNDARY PROBLEMS FOR DIFFERENTIAL 

EQUATIONS* 

BY 

R. G. D. RICHARDSON 

Introduction 

One of the most important problems of mathematical physics is the study 
of boundary problems for differential equations of the second order in one 
or more variables. These equations contain a parameter which must be so 
determined that solutions exist satisfying prescribed linear homogeneous 
relations between the values of the function and its derivative on the boundary. 
The problem in one dimension has been studied in great detail, f As typical 
of the problems in one independent variable we have the equation 

(1) v xx + [r(x) + \k(x)]v(x) = 0, t>(0) =»(1) = 0. 

By means of his theory of integral equations, HilbertJ extended the theory 
of boundary problems in one dimension to the corresponding problem in 
partial differential equations of the elliptic type, of which 

(2) u xx + u vy + (r(x, y) + \k(x, y))u = 0, u = on boundary, 

is typical. By hypothesis, for both equations (1) and (2) the function k has 
one sign (orthogonal case) or k has both signs and r == (polar case). It has 
been shown that there exist an infinite set of real parameter values X* for 

* Presented to the Society, September, 1914. 

t For existing methods and literature of the subject see B6cher, Encyklopaedie der mathe- 
matischen Wissenschaften, IIA7a; Boundary problems in one dimension, Proceedings of 
the International Congress of Mathematicians (Cambridge, 1912), and 
Lecons sur les mHhodes de Sturm (Paris, 1917); Lichtenstein, Zur Analysis der unendlich- 
vielen Variabeln, Rendiconti del Circolo Matematico di Palermo, 
vol. 38, p. 113. 

t Grundzuge einer allgemeinen Theorie der linearen Integralgleichungen, p. 58 (Teubner, 
1912). This book is a reprint of a series of articles in the Gottinger Nachrichten. 

Refinements concerning the nature of the field R , of the coefficients of the equation and 
of the functions / ( x , y) to be expanded in terms of the set of solutions constitute the chief 
progress made in this field since Hilbert published his work. See the literature cited by 
Lichtenstein, Randwertaufgaben der Theorie der linearen partiellen Differentialgleichungen 
zweiter Ordnung vom elliptischen Typus, Journal fur Mathematik, vol. 142, (1912) 
p.l. 
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which solutions of (1) and (2) exist satisfying prescribed boundary conditions. 
In the orthogonal case this set has a limiting point at positive infinity only 
and in the polar at both positive and negative infinity. 

The boundary problem for partial differential equations of the hyperbolic 
type so far discussed has been of an entirely different nature. The existence 
of a unique solution has been shown* in case the value of the function is given 
on one of two intersecting curves and the value of the normal derivative on 
the other, neither curve being cut twice by a characteristic. The elliptic 
and hyperbolic types of equations are essentially different and are char- 
acterized by Du Bois Reymond as " himmelweit verschieden." On the 
other hand, in certain special cases at least, a boundary problem for the latter, 
analogous to that for the elliptic case can be solved, and functions / ( x , y) 
may be developed in terms of the solutions corresponding to the various 
characteristic numbers. A simple example is given by the solutions of the 
equations 
[3(a)] u xx + u yy + \u = 0; 

[3(6)] u xx — Uyy + \u = 0. 

If we demand that solutions of these equations shall vanish on the boundary 
of a square whose opposite vertices are (0, 0), (ir , x) we find that solutions 
sin mx smny ( m , n positive integers) of one equation are solutions of the 
other. For the elliptic equation the characteristic numbers must be of the 
form m? + n 2 , for the hyperbolic, m 2 — n 2 . In the latter case the set has 
limiting points at and at ± w while in the former this occurs at + <» only. 
If the field is taken with sides of incommensurable length there may be an 
infinite number of limiting points for the set of characteristic numbers of 

[3(6)].t 

Including [3(a)] and [3(6)] as special cases there is a class of partial differ- 
ential equations for which known methods for ordinary equations suffice. 
For, if in the unit square we consider the equation 



m ^ ) k{ v{x)d ^- s{x) Ty[ q{y) Ty] 



(4) 

+ [r(x)m(y) — s{x)l{y)\u + \s(x)m(y)u = 

for which solutions are to vanish on the boundary, we know that the functions 
u(x,y) = v{x)w{y) and the corresponding parameters X rai „ 2 = m»i ~ v n% 

* For the literature of the subject and the most general theorem see Mason, Selected 'Topics 
in the Theory of Boundary Value Problems of Differential Equations, The New Haven Mathe- 
matical Colloquium (Yale University Press, 1910). 

t For example, if the opposite vertices are (0, 0) , ( w, ■*/ \/2) , the characteristic numbers 
are of the form m 2 — 2re 2 . Among the integers which occur among this set, each an infinite 
number of times, are all those of the form H 2 — 2 , where H is an integer. 
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will fulfil the conditions, v(x) and w (y) being solutions of the equations each 
in one variable 



d V dv(x)~\ 

;[i(y)— dy ~\- 



[ r ( x ) + ns(x)]v (x) = , 

p > 0, s > 0, »(0) = »(1) = 0, 
[l(y) + vm{y)}w{y) = 0, 



d_ 
dy\ 

q > or q < 0, m > 0, w (0) = w(l) = 

and /x„ i; j-nj the parameter values corresponding to solutions which vanish 
fti — 1 , ra 2 — 1 times respectively. This is immediately evident on multi- 
plying the first equation by m(y)iv (y) and the second by — s(x)v(x) and 
adding. If q > the equation is hyperbolic; if q < 0, elliptic. 

None of the methods at present in use lends itself to the discussion of the 
new boundary problem of the hyperbolic equation. The objects of this memoir 
are (a) to develop a new method by which the well-knotvn results in regard to 
boundary problems for the ordinary and the elliptic partial equations can be 
deduced and (b) to attempt an application of this method to an investigation of 
the facts concerning the new boundary problem for the hyperbolic equation. 

The new method consists in dividing the region into equal rectangles by a 
lattice and studying the difference equations approximating to the differential 
equation at the lattice points. The size of the lattice division is decreased 
indefinitely and it is shown (§ 3) in all cases (including the hyperbolic) that 
each set of corresponding parameter values X has a finite limit. The set of 
numbers which are limits of the characteristic values of the difference has, 
both in the elliptic and hyperbolic cases, a limit point at infinity. In the case 
of the elliptic partial differential equation it is shown (§ 5) that the corresponding 
functions of the lattice approach uniformly a function u(x, y) which is con- 
tinuous and that this limiting function has continuous derivatives of as high 
order as those possessed by the coefficients of the equation. The totality of such 
functions furnishes the characteristic solutions of the equation. 

The results obtained in the hyperbolic case (§ 6) suggest that with the 
possible exception of special cases, the boundary problem similar to that for 
the elliptic case may be possible of solution. However, the reasoning is not 
complete and the problem is left in large measure unsolved. 

Although the problem in one dimension is not treated by the new method, 
its discussion is naturally simpler than that of the elliptic equation which is 
discussed in full. All the ordinary theorems both for the orthogonal and polar 
cases of (1) can be derived in simple fashion by the new processes. The method 
can also be extended (a) to problems with more general homogeneous boundary 
conditions and (6) to problems in three or more dimensions. 

Beside a completion of the discussion of the hyperbolic case, other new 
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and interesting problems of various kinds are suggested by the results obtained. 
No problem in quadratic forms or in linear equations in an infinite number of 
variables in which the set of characteristic numbers has limit points both at 
infinity and in the finite region has as yet been solved. We can however 
set up by a method analogous to that used by Lichtenstein* for ordinary 
differential equations, problems of this nature equivalent to the hyperbolic 
case of (4) for which we know solutions exist. The hyperbolic differential 
equation may be regarded as the Lagrange equation of an isoperimetric 
problem in the calculus of variations in which there are an infinite number of 
integral sub-conditions (see footnote, p. 496). Such a problem has never been 
studied. The expansion of linear forms in terms of the solutions of the 
approximation equations (§ 1, III-IV) suggests an expansion of functions in 
terms of the characteristic functions of our equation. The author hopes 
later to report progress along some of these lines. 

1. The linear algebraic equations of approximation 

Let us approximate by algebraic equations to the most general linear self- 
adjoint homogeneous differential equation of the second order in two inde- 
pendent variables which may be written 

du(x, y) 



d f du{x,y)~\ d I" . 



(5) ftcf '"' dx J dyl^ ""' dy J 

+ [r(x, y) +~Kk(x,y)]u(x,y) = 0, 

where X is a parameter. The notation will be simplest if the region to be 
considered is a unit square with opposite vertices (0,0), (1,1) and there 
are m/x equal rectangles formed by m — 1 lines parallel to the y-axis and ix — 1 
lines parallel to the a;-axis, the values of the functions p , q , r , lc at one of the 
lattice points with coordinates (i/m, j/p) (i = 0, 1, • • • , m; j = 0, 1, 
• • • , m) being denoted by p.-y, la, i"ij, hj, and the derivatives u x , u y being 
replaced by m(«,7 — «<-iy), m( m ;/ — ua-i) respectively. In this way we 
avoid the complexity of notation that arises in the general region for lattice 
points near the • boundary.! For the algebraic equations it will be assumed 

* Loc. cit. 

t More generally the plane may be divided into rectangles by lines parallel to the axes, 
the lattice points denoted by (i,j), the derivatives du / dx at these points replaced by 
( un — Ui-u ) 1 1% where U is the length of the interval between i — 1 , j and i , j , the derivative 
a 2 u I dx 2 replaced by 

2_ / Uj+lj — U,jj _ Uyj — Uj-lj \ 

h+i + U \ k+i h J 

etc. For points near the boundary the lengths of the intervals cannot be made uniform and the 
equations set up for these points must be of irregular form. While in the following discussion 
details are given only for the case that the region is a square it should be borne in mind that 
this restriction is due to a desire to retain simplicity in the notation and that the results hold 
for more general regions. 
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that p, q , r , and k are continuous functions* of x and y and (to insure the para- 
meter values being real) that k is positive throughout the region. For boundary 
conditions we postulate^ that ua vanish on the contour, 

(6) Wo/ = u mj = u i0 = u ifL = 

(i =0,1,2, ■■■,m;j = 0,1,2, ■■-,?). 

There are then (m — 1) (m — 1) homogeneous equations in as many un- 
knowns of the type 

L (ua) + \hj Uij = m 2 [ p i+1 j{u i+lj — Uij) — pij(Uij — «*_iy) ] 

{A) - m 2 [ qij+i (Uij+i - u^) — qn{Uij - Uij-!)] 

+ (r t j + \hj)uij = 

(«' = 1 , • • • , m - 1; j = 1 , • • • , m - 1 ) . 

If we hold fixed the lattice point in question and indefinitely increase the 
number of divisions (in this paper we shall consider in general the set obtained 
by repeatedly doubling m and fx ) the limit of (A) gives the equation (5) for 
this point. 

Intimately connected with equations (A) are quadratic forms 

Q 1 - ^7, f "^ £ Pi+lj(U i+ !j - Uij) 2 
IflfX ij 

(7) - M 2 Z) Qii+i {Uij+i - u^) 2 - X) rn u% ] , 

ij ij 

't 2 ~ Z - 2-1 kij Uij , 

in variables «,-y . From the hypothesis on k ( a; , ?/ ) it follows that Q2 is definite. 
In seeking the minimum of Qi subject to the restrictions (6) and Q2 = 1 , one 
is led on differentiation with regard to each of the variables u t j to the system 
of equations (A). Another relation under the conditions (6) is given by the 
formula 

(8) Qi = - — X) L ( u a ) u u = x & , 

mfx ij 

which is established by multiplying the equations (A) by M,-y, adding, and 
dividing by m\x. . 

In order that there be solutions different from zero of the (m — l)(ju — 1) 
linear equations (^4), X must be a root of the determinant formed from the 

* Eventually the important cases are those where p and q have each one sign. When they 
have opposite signs (p>0,2<0)we shall call the difference equations elliptic and when 
the same sign ( p > , q > ) hyperbolic. 

f More generally it might be assumed that on the boundary h>j «<,- + Uj A„ («;, ) = where 
A„ ( ua ) is a first difference expression corresponding to the normal derivative, and ha , kj 
are sets of constants. The discussion in that case would be similar. 
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coefficients of that set of equations. This symmetric determinant has 
( m — 1 ) ( ;u — 1) roots for X , and from our hypothesis it follows that these are 
all real.* For, since X&*y occurs only as part of each term in the main diagonal, 
we may multiply the row and column in which it stands by 1/ ■\k^j without 
altering the symmetry of the determinant and thus reduce it to the form of the 
secular equation. f Corresponding to each of these (m — l)(/x — 1) values 
of X , called characteristic values, there is a solution u t j of the linear equations 
(A). These solutions shall be said to be normalized when the w's are multi- 
plied by the constant which makes Q 2 = 1 . The normalized solution corre- 
sponding to the characteristic X w and which consists of a set of constants 
each defined for a point ( i/m , j/p ) of the lattice, shall be called a characteristic 
function and denoted by U ( *J . These sets of solutions Uty and the values 
X ( "' [k = 1, •••, (m — l)(ju — 1)] have various properties, some of which 
we shall now develop. 

I. Characteristic functions U ( 'j , Uty corresponding to two values of X are 
orthogonal, that is 

To prove this consider a set of equations (A') similar to the set {A) except 
that pij , qa , Tij , kn , u { j , X are replaced by p\j , q'u , r'a , k'a , u\j , X' respectively. 
Now multiply the equation for u,j in (A) by u\j those for w'.y in (A') by — u i3 -, 
add and sum over ij . There results under condition (6) the fundamental 
relation 

2^i \ r ij \ i^kij r ij X kij)Uij Uij 

lli\X ij 

(10) = - X) (Pa - Pij) ( u ij - u i-u) ( u 'a - u 'i-n) 

M ij 

- — Hilij ~ lij) ( u 'j ~ M *7-l) ( u 'iJ - Uij-i). 

Ill ij 

On setting k tj = k\j, r t j = r\j, pa = pij, Qij = q\j we have for two solutions 
[7$, TJfj of {A) the relation 

(X<"> -X(">) Y.kijZiyuy^O 



* In the problem corresponding to the elliptic differential equation where Qi is definite 
(p > , q < , r <0), the roots X are of course real even if the fc's have both signs. But if 
neither of the two quadratic forms Qi , Qt is definite, at least a part of the X's and of the solu- 
tions will in general be complex. Similarly if characteristic functions exist for the hyperbolic 
case (p>0,g>0) when k takes on both signs, at least a part of the characteristic numbers 
will in general be complex. The minimum and calculus of variations problems corresponding 
will then have no meaning. 

f Scott and Matthews, Theory of Determinants (2d ed.), p. 147. 
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and when w 4= A (v) this gives the desired orthogonal property. If X is 
multiple valued (X (<) =A (r) ), the solutions U ( *} , Uty may by the usual 
methods be so chosen that they are orthogonal to one another as well as to 
the others. The set of (m — 1) (;* — 1) solutions thus chosen as orthogonal 
and normalized is called the fundamental set. 

II. The fundamental set U'f] [k = 1,2, • • • , (m — l)(ju — 1)] comprise 
(m — l)(p — 1), linearly independent sets of (m — l)(ju — 1) constants. 
For, if there exist constants c K such that }£« c K U { "] = for all ij , on multi- 
plication of each fundamental set by & i; - Vf] and summing over ij and then 
over k we get by using the orthogonal properties the equation c p = . This 
equation holds for all values p = 1,2, • • • , (m — 1 ) (ju — 1 ) . 

III. Any set of (m — l)(jt — 1) constants </>;y defined for the lattice 
points is expressible in terms of the sets of constants U ( $ as follows: 



0,7 = Zc U?j ( c K = J- s # fc D« \ . 

* \ m/* v y '■> ,J / 



For, denoting by Ra the remainder <£;,- — £« c,, £7,7 and multiplying by 
kijU^J and summing, we get 12a hi Ra U'-'j = for all k. Since the 
(m — 1)(ju— 1) sets of constants U ( $ are linearly independent by II, the 
(m — l)(ju — 1) multipliers k a it,-y must all be zero. And since by hypothesis 
&*,• > , it follows that Rij is zero. 

IV. The arbitrary set of constants 4>a is n °t generally expressible linearly 
in terms of v of the £7»/s [f < (m — 1) (ju — 1)] which we denote by UfJ. 
However, if we define the coefficients c- to be 

the best approximation to <f>u is given by the formula £j c- Uf) , if by best 
approximation is understood that which furnishes 52»7&»7 [<£»i — ^c- UfJ] 2 
with a minimum. 

For, in order to obtain a minimum for this last expression we equate to 
zero the derivatives with regard to c s and obtain the relations 

T l k ij (<t>ij-i:c,Uf))U^ = 

for all k. From the orthogonal and normal properties it follows that this 
may be written 

C; = X) 0.7 hi U?, 

and on taking the second derivatives these relations are found to be sufficient 
to establish the theorem. 
It follows from III that 4>a — £; c- U'fj = £ s c- U^ where k ranges over 

Trans. Am. Math. Soo. 38 
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the remaining (m — 1) (/* — 1) — v oi the k's and c- is defined in a manner 
analogous to c- . Hence 

And in particular if k is a null set, 

~ Zj *»7 #y = ^ c « • 

V. ^4?i# characteristic number can be considered both as the minimum and the 
maximum for certain problems. If the (m — I) (p — 1) numbers are arranged 
in a system Xi = X2 = • • • X»_i = X» = X s +i = • • • = \ m -i)( ll ,—i) , then X s is 
the minimum* and maximum respectively for those values of Uij which make 
§2 = 1 and satisfy the relation (6) and the relations 

[11(a)] £fcyl7$Mtf-0 («-l, ••.,.- 1); 

[11(b)] Y,kijU ( $Uij = («-«+i, ...,(m-i)(M-l)) 

respectively. 

Since proofs for all cases are similar, we shall content ourselves with a dis- 
cussion of the last case. From III and relation [11(6)] it follows that 

From (12) and (A) we have 

£ w (««) = t,K^k i] -U ( t) 
«=i 

and hence from (12) and (9) 

Qi = -J^L (uij)Uij = Xi c\ + X 2 c| + • • • + X s cf . 

ij 

Then Qi — X 8 = (Xi — X s )c? + • • • + (X„_i — X 8 )c?_i which is zero or nega- 

* As an extension of this result it seems probable that any characteristic function u ( x , y ) 
of the hyperbolic differential equation (5) ( p > , q > ) may be regarded as the solution 

of the problem of finding a minimum and maximum respectively of the integral J (pul 
— qu 2 y — ru* ) dR for those values of u which vanish on the boundary, satisfy J ku* dR = 1 
and an infinite number of orthogonal sub-conditions of the form I kUiitdR — 0, where i 

ranges over the whole set of solutions corresponding to values of X less and greater respectively 
than the one considered. The minimum or maximum value would be the corresponding 
characteristic number X. The minimum property can be easily established for the elliptic 
case ( p > , q < ) , since here the number of sub-conditions is finite; it is probable that 
the maximum property holds also. 



(12) 


Uij 




and hence that Q 2 = 


= c\ + c\ + ■ • 


•+c° = l. 
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tive, and zero when ci = c 2 = • • • = c s _i = . That is, Qi is a maximum 
for U ( ij and takes on value X„ in this case. 

VI. The solution Vf} of V furnishes the maximum and minimum 1/X 8 
respectively for Q 2 when those sets of constants «,,- are considered which 
satisfy the relation Qi = 1, the relations (6) and the relations [11(a)], [11(&)I 
respectively. 

2. Variation properties of the characteristic numbers of the 
algebraic equations 

The functions p(x,y), q(x,y), r(x,y), k{x,y) may be regarded as 
analytic functions of a parameter t in an interval -lSigl. If in the 
discussion of I, § 1 we allow p' , q' , r' , k' to approach p , q, r , k , the param- 
eter X' approaches X and the formula (10) becomes 

^i_Vt „ 2 - _-i.V^'<,/ z _ ^ -r^dkjj . 
dt m M V ' ~ V« dt Uij mfi^f dt Ui} 

From this formula one can determine the method of variation of the param- 
eter X with the variation of one or more of the functions p , q, r , and k . 

I. In the set of equations (^t) all of the characteristic values X are increased 
if one or more of the p's are increased or one or more of the q's decreased and 
all are decreased if one or more of the p's are decreased or one or more of 
the g's increased. 

II. In the set of equations (A) all of the values |X| are decreased if one or 
more of the k's are increased. If X = , the new characteristic value is also 
zero. 

III. In the set of equations (^4) all of the values X are decreased if one or more 
of the r's are increased. 

It may be observed that since 

~~ 2-i *»7 u ij — Qi = 1 > 
mfj, ij 

the change of X is an infinitesimal of the same order as the mean change of 
the r's and k's. And since neither 

— T, m 2 ( u^ - Ui- Xj ) 2 nor ~ , E M 2 ( Uu - «**-i ) 2 
mfi y m/j. ^ 

can be zero when Q 2 = 1 , the change of X is of at least as high order as the 
mean changes of the p's or g's. 
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3. FlNITENESS OF CHARACTERISTIC NUMBERS 

The comparison theorems of the preceding section may be used to determine 
intervals within which the values of X for all m, p must lie. In the unit 
square P with opposite vertices (0,0), (1,1), constants Pi, P 2 , Qi, Q2, 
Pi , P 2 , may be chosen so that 

P lS |>P 2 , Q lS =|<Q 2 , ftglsft. 

From I— III of § 2 it follows that the characteristic values of (A) are inter- 
mediate between those of 

m 2 Pi(u i+ ij — 2un + Ui-ij) — m 2 Qi (Uii+i — Zuij + «,7_i) 
{14) 

+ (Ri + \)u ij = 0, 

m 2 P 2 (Mi+iy — 2M ti + Ut-ij) — n 2 Q 2 {u ii+ i — 2u i} - + u i} --i) 
(15) + (P 2 + X)w t7 = 0. 

That the characteristic numbers of (14) and (15) are the differences between 
the characteristic numbers for the sets of equations 

m 2 Pi(« i+ i - 2v t + vt-i) + (Pi + vi)v { = 0, 

(16) 

ix 2 Qi(w j+ i - 2wj + Wj-i) + ViWj = 0; 

m 2 P 2 (« i+ i - 2vi + Vi-i) + (P 2 + v 2 )v { = 0, 
(17) 

^ Q 2 (Wy+i - 2Wj + Wj-i) + Vi Wj = 0, 

respectively, where v t , Wj are each sets of variables in one dimension and 
v\, v{, vi, v' 2 are parameters, may be seen on multiplying the first of each 
set by Wj and the other by — v t , adding and setting u^ = v t Wj. The X of 
(14) is vi — v\ and of (15) is v 2 — v' 2 . 

It is easily shown that characteristic numbers of (16) are approximations 
for the corresponding differential equations 

(18) Piv xx (x) + (Pi + vi)v(x) =0, „(0) =*>(1) =0, 

(19) Qiw yy (y) + v[w(y) = 0, w(0) = w(l) = 0, 

each in one variable and with constant coefficients. These have solutions of 
the form v = sin «i ttx , w = sin n 2 iry ( n\ , n 2 = 1 , 2 , • • • ) , which oscillate 
ni — 1 and n 2 — 1 times respectively; the corresponding parameter values are 
vi = (niwyPi - Pi and v[ = (n 2 7r) 2 Qi. If we multiply (18) and (19) 
by w(y) and — v(x) respectively, add and set u(x, y) = v(x)w(y) we 
obtain the differential equation 

(20) {Pxu x ) x - (Qiu y ) v + [Pi + (vi-v[)]u = 
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corresponding to (14) and for which explicit solutions sin n\ irx sin n 2 wy 
may be written down. The parameter value X of (20) is then the difference 
v\ — v[ = (wi 7r) 2 Pi — jBi — (w 2 ir) 2 Qi between the parameter values for 
the equations (18), (19), and is characterized by the number of oscillations 
n\ — 1 , n 2 — 1 of their solutions. Using the notation X^, they may be 
ordered as follows: Xu, X21, X12, X31, X22, X13, ••• . Corresponding to the 
first we have a solution sin irx sin wy of (20) which does not oscillate, to the 
second sin 2-kx sin wy which oscillates once as a function of x , etc. 
The solutions of the equation 

(21) k[P 2 u xx -Q 2 u uy + (R 2 + \)u] = 

corresponding to (15) are again of the form sin n\ irx sin n 2 wy but now the 
characteristic number X„ in2 is ( n\ ir ) 2 P 2 — P2 — ( n 2 w ) 2 Q 2 . 

For any given n\ , n 2 and an arbitrarily chosen e > , we may choose m , /I 
such that the characteristic numbers of (16) and (17) differ from those of 
(18) and (19) by less than e/2 when m S m , /* =S fi . Hence the characteristics 
X^nj of (14) and (15) differ from those of (20) and (21) respectively by less 
than e . There is then a parameter value of (A) which lies between 

7r 2 (n?Pi -nlQ 1 ) - Pi + e and w 2 (nl P 2 - n\ Q 2 ) - R 2 - e. 

The parameter value \ ni „ 2 of the general equation (A) shall be considered as that 
which is obtained by variation when the coefficients of (14). or (15) are varied 
continuously in order to assume those of (A). 

The uniqueness of this determination of the X's may be shown by the 
following argument. Let p , q , r be chosen as monotone functions of t in the 
interval — 1 Si t Si 1 so that with increasing t the constant coefficients P ( t ) , 
Q(t), R(t) vary monotonely from P 2 = P(l), Q 2 = Q(l), R 2 = R(l) 
to Pi = P(- 1), Q t « Q(- 1), Pi = R(- 1): the parameter \ ni „ 2 (t) 
will increase from the value X^,, ( t 2 ) = ir 2 ( n\ P 2 — n\ Q 2 ) — P2 to the value 
XnjMji ( h ) = ir 2 ( n\ P\ — n\ Q\ ) — Pi . If the coefficients remained constants 
throughout, for each value of t the solution would be the product of a function 
of x alone oscillating n\ — 1 times by a function of y alone oscillating n 2 — 1 
times. 

The value of the parameter will decrease if the coefficients vary monotonely 
in any fashion from kP 2 , kQ 2 , kR 2 to kPi , kQi , kRi even if p(t) , q(t) , 
r(t) do not remain constant throughout the region for intermediate values 
of the parameter t. If the values of the coefficients differ at most only in- 
finitesimally from constants for any value of t , the value of the corresponding 
X, in2 will differ by at most an infinitesimal from the value X ni „ 2 (t) . By the 
introduction of another parameter 5 to take care of the variation of the coef- 
ficients from the original values P (t) , Q (t) , R (t) , we can have X„ lK!! ( t , s ) 
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a decreasing function of t for all values of s and a continuous function of s for 
all values of t. But the equation (A) may be regarded as obtained in this 
manner by a monotone variation of the coefficients p{x, y, t, s) , q(x, y , 
t, s) , r(x,y,t,s). Since each of the coefficients p ( 1 , s ) , q (1 , s) , 
r(l,s),p( — 1 , s) , q(— 1, s) , r(— 1 , s) is a constant for all values 
of s and since Xn in2 ( 1 , s ) , X ni7l2 ( — 1 , s ) are uniquely determined by n\ and n 2 
it follows that X ni „ 2 ( t , s ) is completely determined in this manner. 

Theorem I. For each value of m, tx there is a characteristic number \i in2 
of the equation (A) determined by the oscillation numbers n\ , n% of the solution 
of (14) or (15). For fixed tii, n 2 the totality of Xn in2 lies within a finite interval, 
upper and lower bounds for which may be readily determined. 

Each set of X„ in2 's will have at least one limit number, and we thus can 
pick out a denumerable set of limit values for the characteristic numbers. 
The bounds of each characteristic number so obtained may be readily calcu- 
lated and definite statements made in regard to the limit points of the set 
for the elliptic and hyperbolic cases. For the former (p > 0, q <0), there 
is a limit point at positive infinity only; for the latter (p > , q > 0) , there are 
limit points at both positive and negative infinity beside possibilities of more in 
the finite region. 

4. Some auxiliary theorems 

For subsequent use we shall now derive some miscellaneous auxiliary 
theorems relative to maxima and minima of quadratic forms and integrals 
and to the limits of linear forms. 

Theorem II. The value of the integral of x(du/dx) 2 dx along the straight 
line joining the points Pi ( , M ) , P 2 ( 1 , N ) ( N being less than M ) is less than 
that along any other curve lying above the line and joining these points. 



P l (0, M ) 




P 2 (1, N) 




Fig. 1 



To prove this we shall first show that the integral taken along two straight 
lines Qi Q 2 , Q2 Q3 is greater than along Qi Q s (Fig. 1). Without loss of gen- 
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erality it may be assumed that the ordinate of Q 3 is . We wish to prove 
that 

£ x {^~) dx < £ x {f 2 ^ ^ + £ x {-^) dx 



or 

at + a, < ^+^ 1(0i _ A) , <* + «» 2 

as — <X\ «2 — &1 «3 — «2 

which may be written 

01 ( « 3 — « 2 ) , 02 ( «8 — ai ) . , 

5-7 C OLi + -5-7 r <X 2 > «2 + OCl . 

jS 2 ( «3 — «1 ) Pi ( «3 — «2 ) 

But if a 2 > en > and < / < 1 we know that fa\ + a 2 // > a 2 + «i . 
Hence our inequality is satisfied when pV/3 2 < (0:3 — oti)/(a 3 — a 2 ) , which 
certainly holds when Q 2 lies above the line Qi Q 3 . 

Since the integrand is never negative the integral along the dotted line 
Ri R 2 whose slope is zero is less than that along the curved line Ri R 2 . In 
the discussion it is then necessary to use only monotone comparison curves. 
But such a function can be approximated as close as we wish by a monotone 
function composed of broken lines. By repeated use of the argument made 
above, our theorem is established. 

Analogous to the division of the plane into rectangles as outlined in § 1, 
one can consider an interval 0, 1 in one dimension as being divided into m 
parts and at the points of division i/m (i = 0,1, ■ • • ,m) variables »; as 
defined. Within the intervals the function v may be considered as made up 
of straight lines joining the points ( i/m , Vi ) . For convenience sake the 
intervals will be taken equal in length. 

Theobem III. 7/ within a subinterval a/m, {a + a')/m the oscillation of 
the Vi's is M, then m Y.1=Z' («i - »i-i) 2 > mM 2 /a' . 

For, we have m Uli?' (Vi — »t-i) 2 = m/a' [ Y,"="'\vi — «<-i| ] 2 and since 
£?=''!«,■ — Vi-i I is at least equal to the oscillation, the theorem is established. 

Cobollaby. When m 2;(?>; — fli-i) 2 remains limited with increasing m, 
it is impossible that the v's have a limiting function V (x) which is discontinuous. 

For, at a point of discontinuity of V let M be the oscillation. Given e > , 
it is possible to find an m such that in a subinterval of length cr'/m < e the 
oscillation of the approximating function is greater than M/2. Then 
™ £i(»t — »i-i) 2 > M 2 /4e which can be made great at pleasure. This 
being a contradiction, V (x) must be continuous. 

Returning now to the discussion of the two dimensional difference equa- 
tions (^4), we shall for the sake of convenience set m = ix and denote by A m 
the maximum of the normalized solution for a given X„ ini! and a given m. 
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Let A m be taken at a point (<r, r) [u (<r, t) = ^4 m ] , and denote by (c + <r', 
t + t') the point furthest from {a , t) such that there is a chain of points 
connecting them for which w,y > A m /2 . 
Theobem IV. If A m is the maximum* of those values utjfor which 

(22) iZ4 = l, 

we have the formula 



mf 



1^ ., .. . .. .._ At'' 2 



O" 



(23) — 2 2 m2 ( M w - M *-W ) 2 + ™ 2 ( «•>• - «.j-i ) 2 = 32 2 m 2 • 

In proving this it may without loss of generality be assumed f that a' i= r' 
i£ . From (22) it follows that the total number of points where Uij > A m /4, 
is less than l%m 2 /A 2 m . Let us select a rectangle with opposite vertices {a , r) 
and (<t + a' , t + 32m 2 /A 2 m <r') having double this number of points in it. 
On at least half the lines parallel to the y-axis in the rectangle, there is then a 
point where u^ s= A m /k and on every such line a point where u { j is A m /2 . 
From Theorem III we have 

J_V *r _ v ^la'^(A m \Ala' Aj<r' 2 

m 2^ m i u a u a-i) = 2 my \ 4 ) 32m ~32 2 m 2 ' 

from which (23) follows at once. 
Were then 

,. Ala* 

hm = oo , 

m=oo ^ 

the value of (23) would increase without limit as A m increases. We can 
investigate the conditions under which 

A 2 «' 

lim^* * oo. 

m—.w m 

Since in this case a'/vi is an infinitesimal of order at least as great as 1/A 2 m , 
and since the same result would be obtained if we used fractional parts of A m 
other than A m /2 , A m /4 , it follows that the point where u t j increases indefinitely 
is isolated. 



* A similar theorem can be obtained for the minimum by replacing m,-,- by — i*,,- . 
f Were there not at least one such point on one of the four directions from ( a , t ) such 
that un > A m 1 4 we should have two particular terms, viz., 

—Am?(u ar — M a _i T ) 2 + m 2 (Mo-T — «<tt-i) 2 ] 
mr 

of the sum (23) greater than \A 2 m but since from (22) A\ % < m? , it follows that this is greater 
than A^ / 2m 2 ; this relation satisfies all the purposes of the theorem. 



1917] BOUNDARY PROBLEMS FOR DIFFERENTIAL EQUATIONS 503 

Theorem V. If at a point P a function u defined over a unit square R takes 
on a maximum M and if 






u 2 dR = l, 
then 

r {ul + ul)dR ^M*o* M / 32 2 , 



I< 



where u M is the. radius of the greatest circle such that there is a point Q on the cir- 
cumference which may be joined to P by a curve along which u S M/2 . 

The proof of this theorem is analogous to that of Theorem IV. 

Let us now deduce some properties of solutions of the elliptic partial dif- 
ferential equation which for the sake of simplicity may be taken in the form 

(24) u xx + u vv + (r + \k)u = . 

If this is multiplied by u and integrated over a field S bounded by a curve C 
there results the formula 

I ( r + \k ) u 2 dS = I ( ul + u\ ) dS + I u-^-ds , 

where du/dn is the normal derivative directed inward; and hence 

(25) max(r + \k)fu 2 dS^f(ul + u 2 v )dS + f u^~ds . 

Using the notation of Theorem V, let K denote a circular conical surface 
having for its axis the z-axis and cutting the plane z = M/2 in a circle of 
radius <r M . The surface u(x, y) lies outside the cone at the vertex and 
inside or on the surface at the section z = M/2. Denoting by C\ the pro- 
jection on the xy-plane of the intersection of K and the surface u(x, y) and by S 
that part of the region R in the xy-plane bounded by d , then 

(26) f(ul + ul)dS >t^- fu 2 dS. 
Js **GmJs 

For, on passing to polar coordinates we have 

(27) f(u 2 x + ul)dS= ftul+^pdpdd^ fulpdpdd. 

'JS 'JS \ P / Js 

From Theorem II we see that this last integral is greater than the value 
computed as if u were on the conical surface, in which case u^ would be 
- M/2<x M . Hence from (27) 

(ut + u 2 y )ds>^J P 2 0l {e)dd, 
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where p Cl is the value of p on the curve C\ . But since u Si M it follows that 

(29) fu 2 dS S i¥ 2 fdS = M 2 f pdpd» = ~ f ' P 2 Ci (d)de. 

Js Js Js " Jl> 

Formula (26) follows now from (28) and (29). 
Theorem VI. The value of the integral 



I (ul + ul)dR 
Jb 



for a normalized solution u(x, y) of the elliptic differential equation (24) having 
a maximum M is greater than 

M 2 
64 2 max (r + \k)' 

For, on the curve Ci the derivative du/dn is positive and hence the line 
integral in (25) is positive. Using formulas (25) and (26) we get 



(30) 



f(u 2 x + ul)dS <max (r + M) f u 2 dS 
Js Js 



< 4<tm max (r + \k) | (ul + u\ ) dS , 
Js 



and hence 



4<r2jf > max(r + X^) 



From this formula and Theorem V we have 

M 4 A M 4 



r m a m 

(31) J B ( u ° + ^ dR> ^>w^^(7+W)- 

Corollary I. The value of the integral j ( ul + u 2 y ) dR is greater than 

Jb 



M 4 



64 2 max (r + M) f u 2 dR 
Js 

This may be seen by following through the steps of the reasoning when the 
hypothesis that the solutions are normalized is omitted. 

Corollary II. The maximum value of a normalized solution of the equation 
(24) is less than 8 Vmax ( r + X& ) . 

This follows from (31) and (25). 

Let us next consider some properties of the hyperbolic equation and its 
approximating difference equations. For the sake of simplicity of notation 



1917] BOUNDARY PROBLEMS FOR DIFFERENTIAL EQUATIONS 505 

the equations will be taken in the canonical form 

Uxy + (r + \k)u = 
and 

mju {uv — Uij-i — Uj-ij + «i_i.,-_i) + (»"« + \k i:i )uij = 0, 

the solutions being supposed as before to vanish on the boundary of the 
region. The characteristics are now parallel to the axes and by a simple 
summation of the difference equations we get the formula 

(32) [u(A) +u(C) -u(B) -u(D)] = — £(*/ + U i:i )u i} -, 

where A, B, C, D are the vertices of a rectangle with characteristics for 
sides, the summation extending over this rectangle. 
Theorem VII. If for limited Uij the single sum* 

-Yi m2 ( u i} - Wi-lj) 2 

7Tb $ 

is unlimited for some section parallel to the x-axes it is unlimited for all sections 
and the double sum 

— T,m 2 (uij -u^u) 2 
mix ij 

is also infinite. 
For, if we square the special case of the formula (32) 

1 ' 
m(u st — u s - lt ) = m(u, T — u s - lT ) H (r»j + X& SJ )w SJ -, 

where ( s , t) is any lattice point and ( s , r ) a point such that t > r , and then 
divide by m and sum over s, we obtain the relation 

12 m(u it - w,_i<) 2 — Hm(u iT — u^. lT ) 2 

i i 

= ^E/f (r, ft) + E 2/i(r, ft) (u iT - m_ lT ) , 
m i i 

where the totality of the functions 

1 ' 
fi(r, ft) = X) (r ti + \h S j)u si 

M j = T 

is limited, each being the average of a set of limited functions. The first 
sum on the right being the average of the f\ is then limited. If the second 
sum is limited, the sums ]£» m(uu — w*_k) 2 and Si m (u t7 — Wj_i T ) 2 in- 
crease indefinitely together, their difference being finite. If the second sum 
is unlimited, its order of infinity is not more than half that of the sums on the 



' Corresponding to the single integral f u%dx. 
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left. For, we have 

\12fi(r, k)(u iT - m<_, t )| < max |/i|]CK T - M;-i T |, 



Y, \u w - Ui- lr | ^JmJ2 (u ir - Ui- lT ) 2 

i \ i 



Henee the argument holds also for this case. And further if we sum over 
all values of t and multiply by 1/p we are again taking an average and the 
double sum 

— -Y,™ 2 (Uij- Mi-v) 2 
approaches infinity. 

Similar theorems may be readily established for the sums 

X) M ( My — M,-y_i ) 2 and — X M 2 ( «</ - «.-y-i ) 2 • 

In discussing the effect on the solutions of the equations (^4) of increasing 
m , jj. , let us consider the hyperbolic case (which presents the greater difficulty) 
and for simplicity of notation set p i} - = q i} - = 1 and m = p. Adding the 
equations for M»+iy , M»_iy , w,y+i , ua-\ we have 



*2J 



m 2 (u i+ 2j + Ui- 2 j — Uij+2 — Uij-i) + 4 (r,-y + X£,-y) w,-. 

(a) = 4(rij + \kij)uij — (r i+lj + X& i+ iy ) u i+ y — (»v_iy + X&;_iy)w;_iy 

— ( r.-y+i + Xfc»7+i ) Mj-y+i — ( r,-y_i + X&;y_i ) w»y-i , 

while the set of equations for a value of m half the original is 

(/3) m 2 (M i+2 y + Mi_ 2 y - M,-y+2 - «.■,-*) + 4 (r,-y + \hj)un = 0. 

For each value of m the equations (^4) will have a set of characteristic num- 
bers X. Let us choose a X for some m and for larger values of m pick out 
(by the method given in § 3) X's corresponding to this first one, the totality 
of such X's being limited. On the assumption that r and k and their first 
difference quotients m ( ra — r;_iy ) , m ( r<y — r;y_i ) , etc., are limited and 
that the values of w;y and the first difference quotients m(«,-y — w*_iy), 
m ( Uij — Uij-i ) are limited, there will then be an upper bound M for the 
absolute value of these functions and their difference quotients. 

The left-hand members of equations (a), (/3) are identical in form and we 
shall now show that the right-hand side of (a) is in all cases an infinitesimal 
depending on 1/m . For, we note that the expression tij ua — r»+iy M»+iy is in 
absolute value not greater than 

2M 2 



I M»y+i | + | fij || Uij+i — u^ | < ■ 



m 
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We have also 

2M 2 , , 2M 2 



| 'ij Uij *ij — 1 wij — 1 | ^. 



to ' i - - ■■"-•■"-'• m 
2M 2 



m 



and similar relations for the four corresponding terms in k . Hence the right- 
hand side of (a) is less in absolute value than 8M 2 /m [ 1 + upper bound of |X | ] 
and is then an infinitesimal of order at least as high as 1/m. 

As is immediately evident, the same discussion would show that if to were 
multiplied by any power of 2 and the corresponding equations were added 
together in the proper manner, there would result a set of equations each 
differing from (/3) only by terms on the right which taken together form an 
infinitesimal depending on 1/m. Such a state of affairs would lead us to 
suspect that we can for the various values of to pick out corresponding solu- 
tions Ua which are very close to one another and have one and only one limit 
function as to increases indefinitely. 

5. The elliptic equations 

In this section we shall consider the elliptic case of the difference equations 
(^4) (p > , q < 0) . It will be found possible to trace for various values of 
to, n the relations between the corresponding values of X and u^ for these 
equations. When m or n is doubled, the number of equations and the number 
of characteristic values and solutions is approximately doubled. To be more 
exact, it is changed from (m — l)(/x — 1) to (2m — 1) (ju — 1) or 
(m — l)(2jtt — 1). Some of the new characteristic numbers are greater 
and some are less than any of the original set, as one may see by studying 
the minimum problem. There are, however, (m — l)(ju — l)of the new 
set to each one of which there is related one of the old set in that the latter 
is a more or less rough approximation to it and the corresponding members 
of the set of U{/& are similarly related. If m and ix are further increased, one 
may trace the new X's and Ui/s which are related to those already selected. 
And these approach as limit the parameter values and solutions of the elliptic 
partial differential equation 

(33) (pu x ) x ~(qu y )y+ (r + ~\k)u = 0, p > 0, q < 0. 

The details of this discussion we now proceed to give. 

We have already shown in § 3 that for fixed n\ , tii the set of parameter 
values is limited for all to , n . From this it is possible to show that the various 
sets of Ui/s corresponding to the Xn inj 's are limited. With this end in view 
let us prove the following 
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Theorem VIII. If the maximum of the normalized solutions of (A) increases 
without limit with m , /t , the sum 

— £ m 2 {Uij - Mi_iy) 2 + jU 2 (U i} - - Uij-i) 2 

mfi y 

increases without limit. 

To prove this let us make use of the results of Theorem VII and its corollaries 
by passing a surface consisting of portions of planes through the points ( i/m , 
i//A , Ui 3 - ) . Let us agree to pass planes through the sets of three points cor- 
responding to (i,j), (i + l,j), (i,i+l), and others through those 
corresponding to (i + 1, j) , (i,j+l), (i + 1, j + 1) and compare the 
integrals 

fu 2 dR, f (ul + ul)dR 

Jr Jr 

for the total surface found in this way with the sums 

— Z "I > — Z m2 ( u a ~ «*-v ) 2 + M 2 ( u i} - - u^^ Y 

m/A ij TrljX (j 

respectively. 
As a first result let us show that the integral 



f(ul + ul) 
Jr 



(u 2 x + ul)dR 

*JR 

is equivalent to 

(34) — £ m? (Uij - Ui-xj ) 2 + M 2 («« - «.7-i Y • 
m\x ^j 

If we consider a triangle on the surface defined by the points (i, j),(i — 1 , j) , 
(i,j — 1) the slope in the a;-direction is m(Uij — Wj_iy) and in the y-direc- 
tion is ix(uij — Wiy-i) . The integrand being a constant and the area being 
l/2mn, the integral has a value 

(35) 7j— [ m? ( u^ - u^xj ) 2 + m 2 ( u i} - Uij-i ) 2 ] . 

For the triangle determined by (i,j), (i — l,j), (i — 1, j + 1) , the 
integral has a value 

(36) 2^~ t m2 ( u a ~ u i-v ) 2 + M 2 ( Ui-ij+i - Mi_iy ) 2 ] , 

and for the triangle determined by (i,j), (i,j — 1), (i+ l,j — 1), the 
integral has a value 

(37) 2^~ f m2 ( M »+ii ~ u a ) 2 + ^ 2 ( w »7 - M »7-i ) 2 ] • 
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By adding (35), the first part of (36) and the last part of (37) we get a term 
of (34) and similar considerations for the other triangles establish this first 
result. 

As a second result let us show that the sum 

/ j tin 
is greater than the integral 



\b 



dR 



where u is the surface consisting of planes. The integral of u 2 over one of the 
triangles is not greater than 

~max( Uij ) 2 

which is less than 

2^ ( I Ui, | + | ut+ij | + | u ij+1 1 ) 2 . 

Summing over all the triangles we see that this is less than 

V™* 2 

From these two results and Theorem VI, Corollary I, we have the formula 

1 M 4 

— - 23 m * ( u a - u i-u ) 2 + M 2 ( uu - Uij-i f > 



Mii 192=max-(r + X/0^-E< 

where M is the maximum of «,-,■ . Theorem VIII is then established. 
In formula (8) which when written in full is 

— -£ [m?Pij(Uij - Ui-ijf - m 2 qaiuij - w,y-i) 2 ] 

(38) mM « 1 

= ~ 2-j ( r a H~ ™ij ) u ij > 

let us consider normalized values of u { j , ( S»v kij i4j = mu ) , and note that 
since max |r.v|/&»7 is finite and, since for any particular suite of X ni „ 2 on 
which attention is fixed the value of X is limited, the right-hand member is 
limited for all m, p. However, since pij and — qt,- are greater than some 
positive constant, it follows from Theorem VIII that if max Uij increases 
without limit the left-hand member of (38) increases without limit. This 
would establish a contradiction and hence the 

Theorem IX. Those solutions of the elliptic case of equations (A) which 
correspond to a particular set of X„ in!! are limited for all m, p. 
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In proceeding to prove that the first difference quotient m(Uij — M,_ij) 
is limited, we shall confine* ourselves to the case where p and q are constants 
and denote them by x and p respectively. It will naturally be necessary to 
put some restrictions on the first difference quotients of the coefficients of the 
equations. Let us assume that the first derivatives of r and k exist and are 
limited. On differentiation of the equation (24) with regard to x one obtains 

(39) Tru X xx — pu X yy + ( to + \ku) x = 0. 

We note that u = on the whole boundary of the square, u x = and u xx — 
on the sides y = , y = 1 , and u v = and u yy = on the other sides. Hence 
for the special case p = ir , q = p it follows from (24) that u xx = and u yy = 
on the whole contour. If we multiply (39) by u x and integrate over the 
square we note from the resulting relation 



(40) 



f(Tul-pu 2 xl/ )dR= f[(r + \k)u 2 x +(r + \k) x u x u]dR, 
Jr Jr 



that if u x is infinite the order of infinity of the left-hand integrand is higher 
by two than the right, which would lead us to suspect a contradiction. 

With this process for differential equations as a guide let us set up by 
subtraction of the equations for i , j and i — 1 , j in (A) a new set of difference 
equations corresponding to (39). For the case p{x,y) = x, q(x,y) = p, 
after multiplication by m these become 

m 3 X [ (Ui+ij — 2uij + Ui-u) — (Uij - 2ui-ij + u^,-)] 

— mix 2 p [ (u ii+ i - 2iHj + ttiy_i) - (Ui-ij+i — 2«i_iy + Wf-i^i) ] 
(41) 

+ m[(rij - r^u) + \(k i} - la-i,-)] u^ Xj 

+ m ( r,-y + \kij ) ( Uij — Ui-u ) = . 

On multiplication by u i} — w;_iy, summation over i, j, and division by /* 
we get 

X) - ir [ m ( u i+ ij - u tj ) - m ( Un - Ui-y ) f - — p [ m ( u ij+1 - u tj ) 



v M 



m 



(42) - m ( «i_iy - Wf_i.,-_i ) f = - X ( r a + ^u ) ( u u ~ u i-v ) 2 

+ - 13 [ ( r ii - fi-V ) + X ( )fe„ - fei-iy ) ] ( Uij - Ui-ij ) Ui-lj . 
M i} 

Having already shown that the definite quadratic form on the left of (38) is 
limited for all m , ju , it follows that 

my-, 

£_, \Uij Ui—ij } 
M ij 



'■ This change is made here solely to avoid complicated notation. 
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is limited. From this we may deduce that since r.-y + Xfc,-/ is limited, the 
first sum on the right of (42) is limited and further from the formula 

mp. ]T (Uij - Wi-iy) 2 g [ Z) K; - «i_iy| ] 2 , 
or 

— 2 («« - ««-i;) 2 = ~H\ u u - u i-u\ 
that 

is limited. The second sum on the right of (42) being less in absolute value 
than 

max m | ( r iy — r 8 -_iy ) + X ( k i} - — &,_i,- ) | max | ut-y | - 22 1 M »i — M »-i; I > 

must then remain limited for all m, \x. 
The equation (39) may be written 

(irUx)xx ~ (pUx)yy + (f + Xfc ) « x = - ( ^ + X^ ) « 

and regarded as a non-homogeneous equation of the second order in u x . 
All that is essential in the reasoning of Theorems VI, VII, VIII may be applied 
to this equation and the corresponding difference equations when u and u t j 
are limited. Regarding the difference quotient as unit instead of ua itself, 
it follows at once from a theorem analogous to VIII that the left-hand side 
of (42) increases without limit with m ( «,-,- — M»_iy ) . This would lead to 
a contradiction. We have then 

Theorem X. If the coefficients of the elliptic case of the difference equa- 
tions (A) are values taken from functions which have continuous first derivatives, 
the first differences m(uij — w;_iy) , p. ( u a ~ Ua-i ) which correspond to a 
particular set of X„ l7l2 are limited for all m, fx. 

By setting up difference equations by subtraction in (41) or (^1) it may be 
proved in the same way that the three second difference quotients corresponding 
to u xx , u xy , u vv are limited. To indicate briefly the method used in one of 
these cases let us select that corresponding to u xy and confine ourselves to 
derivative notation, remembering that every step may be duplicated by 
difference equations. On differentiating (39) with regard to y we get 

(43) lTU X y XX — DUxyyy + ( Hi + X/Wi ) X y = 

and on multiplication by u xy and integration, the formula* 

, x fl™%*-pulJdR= C(r + Xk)n%dR+ f[(r + U) x u v 

(44) «Ak Jr Jr 

+ (r + \k) y u x + (r + \k) xy u]u X ydR. 

* The vanishing of u xx and u yy on the contour entails the vanishing of «,„ and u xyy also. 
Trans. Am Math. Soo. 33 
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Since the right-hand side of (40) is limited it follows that the left-hand side 
and hence that the integral of u% is limited. The first integral on the right 
of (44) must then be limited and on the assumption that r xv , k xy are limited 
the second is a fortiori limited.* But were u xy itself unlimited the integral 
on the left could be proved unlimited by an extension of Theorems VI, VII, 
VIII, where u xy (and not u ) is regarded as the unit. This would establish a 
contradiction. 

And by an extension of the method, higher difference quotients can be 
shown limited if enough derivatives for r and k are assumed to exist. 

From well-known theorems it follows that since u and its first difference 
quotients remain limited for all m , n , there is at least one limit function u(x,y) , 
which is approached uniformly by the solutions of the difference equations. On 
the assumption that the functions p, q, r, k have limited third derivatives, the 
function u also possesses continuous first and second derivatives. On account 
of this continuity the function u{x,y) satisfying (33) for an everywhere dense 
set of points must throughout the region be a solution of this differential equation. 

From relation (8) it follows that each set of parameter values of the dif- 
ference equations approaches the corresponding parameter value of the 
differential equation uniformly. The X's are roots of an algebraic equation 
of order (m — l)(ju — 1) and being continuous functions of its coefficients 
are continuous functions of the parameter t. Since the Ui/s being roots of 
the linear algebraic equations (A) are continuous functions of the parameter 
t on which the coefficients pa, qa , r^, ku, and X depend and, since they 
approach the function u(x,y) uniformly, it follows that the limit function 
u{x,y) is a continuous function of t . From this we can deduce that there 
is only one limit for the set of Ui/s. For, in the case of equations (16) in one 
dimension with constant coefficients it may be assumed from the theory of 
ordinary differential equations that there is only one limit for each of the sets 
of v's and w's. For the case of the equation (14) in two dimensions with con- 
stant coefficients it follows from this immediately that the %/s have one limit 
function only. In considering as in § 3 the coefficients of (A) as functions 
of t if there is a point t for which there is a single solution U ni n 2 (x,y) of (33) 
while for t + e there are two solutions U„ ini (x, y) , U„ ini (x, y) for every 
positive e , we have as the limit of (9) the orthogonal relation 



I< 



kU ni n,(t + e,x, y)U nini (to + e, x, y)dR = 0, 

which on passing to the limit e = gives 

* The function r + \k with its derivatives and the functions u , u x , u v being limited, the 
only inequality needed is / u 2 xy dB i; J \ u xv \ dR . 
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kUl m (t , x, y)dR = 0. 



L 



This is, however, impossible since the solutions are normalized and the result 
is established. We may now formulate 

Theorem XL There exists in the orthogonal elliptic case (q < 0, k > 0) 
a denumerably infinite set of parameter -values X , with limiting point at infinity 
only, for each of which there is a solution of (33) vanishing on the boundary. 

A similar fundamental theorem may be established for the polar case of 
(33), viz.; when r(i,j)S0 and the region R may be divided into a finite 
number of simple regions in the interior of some of which h is positive and in 
the interior of the remainder of which k is negative. The essential point of 
difference in the treatment of the polar case is in showing that the corresponding 
suite of A ni n 2 's for the various values of m , ix is limited. This may be estab- 
lished by showing that: (a) the absolute value of all parameters X^„, of the 
equations (33) are increased by decreasing the size of the region and retaining 
the condition that u vanish on the boundary; (b) on taking for the smaller 
region a sub-rectangle in which k > , the value of each of the positive X„ in2 's 
is still further increased by replacing k by its minimum (> 0) in the sub- 
rectangle and r by its maximum. This gives an upper bound (which may 
be calculated explicitly) for each positive X ni „ 2 ; (c) on taking for the smaller 
region a sub-rectangle in which k < the value of each of the negative X TC ,n 2 's 
is still further decreased by replacing k by its maximum ( < ) in the sub- 
rectangle and r by its maximum. This gives a lower bound (which may be 
calculated explicitly) for each negative X„ in2 . 

6. The hyperbolic equations 

Preliminary to an investigation of some properties of the hyperbolic differ- 
ence equations let us study some properties of the hyperbolic differential equa- 
tion 

(45) {pu x ) x - {qu v )y + (r + 'Xk)u = (p>0,g>0). 

While an analogous discussion can be given for the more general equation, 
we shall, for the sake of simplicity, confine ourselves to the case where p and q 
are constants. The characteristics will then be sets of straight lines making 
equal angles (7 = arctan ± ^q/p) with the a;-axis. It will be assumed 
that r and k are continuous with as many derivatives as desired and that the 
solutions of (45) are continuous and possess first and second derivatives. 

Theorem XII. If the sides of a parallelogram ABCD are characteristics of 
the differential equation then 

(46) 2Jp~q[u(A) + u(C) - u(B) - u(D)] = f (r + U)udP, 
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the double integral being taken over the interior of the parallelogram. This formula 
holds identically for all such parallelograms and for all solutions of (45) . 




Fig. 2 
For, on integrating the equation (45) over P we get 

I pu x dy — qu y dx = — I ( r + \k ) udP , 

where T is composed of the sides AB , BC , CD, DA. To calculate the line 
integral along AB , we resolve it into two parts, one along AB and the other 
normal to it. We have the formulas 



dx 



i. 



V 



ip + q 

-J — r — - u„ -J- 



ds , dy 



4. 



q 



v + q 



ds , 



Mi ; - — U s \l -— — U„\ " , 

'p + q Vp + q 



4 



j if 

lp + q + Un Vp~+l 1 ; 
and this part of the line integral becomes 

I u s 4pqds = 4pq[u{B) — u{A)}. 
•Jab 



We have further 



I pu x dy — qu v dx — ■\pq [u(B) — u(C)] , 
Jrc 

I pu x dy — quydx — Vpg[t«(D) — u(C) ] , 
Jcd 

J DA 



Formula (46) is then obtained by addition. 

Let us make one or more of the following hypotheses. 
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Hypothesis I. The solutions of (45) are normalized; that is, j ku 2 dR 

Jr 
= constant, where R is a unit square containing the parallelogram. 

If r and X are limited the integral | | ( r + X& ) u | diJjwill be limited. It may 

Jr 

be assumed without loss of generality that this last integral is 1 . 

Hypothesis II. The integral I (ul + u\) dR is limited. 

Jx 

Solutions which satisfy I cannot be infinite at a finite number of points 
in the region, as is readily seen from formula (46). 

Theorem XIII. Solutions of the differential equation (45) which satisfy 
Hypotheses I and II are limited. 

For, let u take on a large value M at the point A . It follows from I that 
there will be points in the neighborhood of A where u < if/2. We may 
distinguish (Fig. 2) two cases: (a) u(B) = u(D) = M/2, while within the 
intervals AB and AID , u > M/2; (b) along at least one of the characteristics, 
u > M/2 . In case (a) we have from (46) and Hypothesis I that u ( C ) < 1 . 
Let 571 and 772 be lengths of AB , AD respectively, and 171 s= 772 . For two' 
corresponding points E and E' we have on applying formula (46) to the 
parallelogram ABE' E , the relation u(E) - u(E') > M/2 - 1. The 
double integral of ul + u\ is independent of a rotation of axes and if the 
a:-axis is taken in the direction AB we have from a theorem of integrals analo- 
gous to Theorem III the formula 

f^-iY 

(ul + ul)dxm ul dx g — • 

Hence 

(K_A 2 

X\o 1 / ~\>T \2 

(ul + ul) dydx g 772 cos 2-y g ( ~ — 1 J cos 27 , 

and since this integral increases without limit with M , a contradiction estab- 
lishes the theorem in this case. 

Case (b) may be readily handled by the methods used in proving The- 
orem IV. 

Instead of allowing u to increase indefinitely we can fix its maximum equal 1 
and have a suite of functions u for which Hypotheses I and II become 

Hypothesis I' lim I u 2 dx = , 

Jr 

Hypothesis II' lhn J ( ul + u\ ) dx - . 

Jr 
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Prom Theorem XII it follows at once that there can be no series of solutions 
of the differential equations having these properties. 

Every solution of the differential equation (45) can be approximated as 
close as we wish in the parallelogram P by a solution of the difference equa- 
tions; by defining a surface made up of portions of planes through the points 
Uij one obtains from equation (A) the formula 

(47) 2^pq[u(A)+u{C) -u(B) -u(D)]=^-J^(u i +U ti )u ii +e mfi , 

TflfJt, tj 

where the sum is taken over the lattice points within or on the parallelogram 
and e m ^ is an infinitesimal depending on m , fi . 

One can, moreover, get a precise formula analogous to (46) for the difference 
equations when p, q are so chosen that the points of the lattice lie on the 
characteristics, that is, if ^p/q is a rational number. Let us derive such a 
formula for the case p = q = 1 . The results in the more general case are 
similar. The difference equations for u { j for the case p = 1 , q = 1 are of 
the form 

Ui+ij + Ui-xj — u ij+ i — Uij^t = — l/mixiuj + \kij)Uij. 

Let us sum these over a rectangle whose vertices A{<r\, n), B(<r 2 , r%) , 
C(<r3, T3) , D(<Ti, Ti) are at lattice points and whose sides make angles of 
ir/4 with the axes. We note that for every point within the rectangle the 
term iiy occurs four times in the sum: once each in a positive sense for those 
equations which enter for points next to the left and next to the right; once 
each in a negative sense for those next above and next below. This is also 
true for those points on the sides of the rectangles which are not vertices. 
For points lying in the next row below AB and to the left, with the exception 
of the end points (<ri — 1, n), (<r 2 , t 2 — 1) , the term occurs twice in the 
sum, once as positive and once as negative. A careful analysis of the remain- 
ing terms gives the formula 

(Vl + M <ri-lTi) + (M<t 3 t 3 + «<r 3 +lT 3 ) - («<r»« + ^Ws-l) 

- ( u <riTi + u airi +x ) = 1/mju X) ( Uj + y<hj ) Uij . 

ij 

Since under Hypothesis II the difference between the values of u at two 
neighboring points is less than a constant 6 , we have 

2[u{A) + u{C) -u(B) -u(D)] = l/mnY,(rij + \kij)uij + ri miJ ., 

( 49 ) * I . A n\ 

In general if VpTtf i s rational we get a formula similar to (48). 
Theorem XIV. When ~*lp/q is a rational number, while 

' / 1 Mij ^ij -*- > 



mn 



V 
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and 

X) ™ 2 ( U U — u i-lj) 2 + M 2 («« — Uij-i) 2 

is limited, it is impossible that max u^ increase indefinitely with m , n . 

This is proved in a manner analogous to Theorem XII. While in the 
case where Vp/g- is irrational precise formulas corresponding to (48) are not 
available, there seems to the author no reason to think that the theorem does 
not hold in that case. This theorem is based on formula (48) which depends 
in no way on boundary conditions but is identically true (when 4p/q is rational) 
for all solutions and all points of the lattice. Further, the solutions of the 
difference equations are continuous functions of the parameter t on which 
the coefficients p, q depend. From the principle of continuity we should 
expect a formula similar to (49) to obtain for all p and q . 

As was seen in § 3, for each m, fi the limited function X„ in2 (t) is monotone 
in the interval —1,1. Each function possesses a derivative and the points 
where — d\ ni njdt == M > form an aggregate whose measure is less than 
c/M , c being a constant independent of M , m, ix. For the infinite set of 
functions Xn,^ corresponding to various values of m, /x there must exist at 
least one limit function and this function will be limited and monotone. It is 
well known that such a function possesses a limited derivative except possibly 
at a set of points of measure zero. 

In formula (13) 

^ _ ^h. J_ V h 2 JL V* ^i 2 A_ V ^ ij 2 

dl = dlm~ix\ ki5 UiJ ~ ~ m M V dt %ii ~ my. 2 * dt UiJ 

+ 4 5 [ m2 1? ( «« - wi -« )2 - ^af ^ - ^- i)2 ] 

we may take (as in § 3) the limited derivatives — dp/dt , dq/dt , dr/dt , dk/dt 
positive and argue that the sum 

— 2~L m2 ( u a — u i-u ) 2 + M 2 ( u a — u a-i ) 2 

is limited for those values of t for which dK/dt is limited. It is then possible 
to enunciate the following 

Theorem XV. For all values of t except possibly for a set of measure zero, 
the sum 

-- £ m2 ( u a - Ui-u) 2 + M 2 ( u n - un-i) 2 
m\x ^ 

is limited. 

If for a value of t for which d\/dt is finite a formula analogous to (48) is 
valid, we can deduce at once that w»y remains limited for all m, (x and that 
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there is a limiting function u(x, y) satisfying the differential equation. If, 
for example, for t = t we have both dk/dt limited and Vp/g rational this would 
be the case. 

From another point of view, any solution of the difference equations is an 
approximation to a solution of the differential equation. With increase of 
m, fj,, a, sequence of limited solutions of the difference equations which at a 
given lattice point are made to take on a value 1 will have at least one limiting 
function which is a solution of the differential equation at that point. But, 
as was seen above (cf. Hypotheses I', II'), for any sequence of solutions of the 
differential equation it is impossible that 

lim I u 2 dx = 0, lim I (vl + vl)dx = 0. 

But this seems to imply that for limited dk/dt, normalized solutions of the 
difference equations like normalized solutions of (45) remain limited. 

It may be proved by a method similar to that used in Theorem VII that 
when the double sum 

remains limited as m, /* increase, the .single sums Z)» m (tiij — M;_iy) 2 for all j 
and 2 3 - /x ( u i3 - — W;y_i ) 2 for all i remain limited. It follows from the corollary 
to Theorem III that the limiting function u(x,y) of the sets of solutions of 
the difference equations (A) is a continuous function for all lines parallel to 
the axes which belong to the lattice for any m, (x and hence is continuous 
throughout the square. 

If we make the hypothesis (which appears a most reasonable one) that a 
formula analogous to (48) can be obtained for any set of hyperbolic difference 
equations (^4), there would be for all values of t (except perhaps for a set of 
measure zero) a bounded limiting continuous function u ( x , y) which would 
be a solution of (45). If there were a value of t for which such a function u 
did not exist, there would be coefficients as near to the original as desired 
such that for these the differential equation would have a solution. 
Brown University 



